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ABSTRACT 

We present a new estimate of foreground emission in the WMAP data, using a 
Markov chain Monte Carlo (MCMC) method. The new technique delivers maps of 
each foreground component for a variety of foreground models, error estimates of the 
uncertainty of each foreground component, and provides an overall goodness-of-fit mea- 
surement. The resulting foreground maps are in broad agreement with those from 
previous techniques used both within the collaboration and by other authors. 

We find that for WMAP data, a simple model with power-law synchrotron, free-free, 
and thermal dust components fits 90% of the sky with a reduced of 1.14. However, the 
model does not work well inside the Galactic plane. The addition of either synchrotron 
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steepening or a modified spinning dust model improves the fit. This component may 
account for up to 14% of the total flux at Ka-band (33 GHz). We find no evidence for 
foreground contamination of the CMB temperature map in the 85% of the sky used for 
cosmological analysis. 

Subject headings: cosmic microwave background — cosmology: observations — diffuse 
radiation — Galaxy: halo — Galaxy: structure — ISM: structure 


1. Introduction 

The Wilkinson Microwave Anisotropy Probe ( WMAP ) produces temperature and linear polar- 
ization radio maps at five frequencies with 1° or better resolution and tightly constrained systematic 
errors. The frequency bands are centered on 22, 33, 41, 61, and 94 GHz; denoted K, Ka, Q, V, and 
W, respectively (see Page et al. 2003 for details). While designed to measure the cosmic microwave 
background (CMB) radiation it also observes the large-scale structure of our Galaxy at angular 
scales and frequencies that are relatively unexplored. Study of our own Galaxy has had a significant 
effect on our understanding of galaxies in general. 

Radio emission from galaxies is generally understood as arising from three effects: “non- 
thermal” synchrotron emission from relativistic electrons spiraling in large-scale magnetic fields, 
“thermal” free-free emission from non-relativistic electron-ion interactions, and emission from vi- 
brational modes of thermal dust grains. At lower radio frequencies the synchrotron emission is 
usually dominant, with flux decreasing at higher frequencies approximately according to a power 
law 14 (f3 « —3). Free-free emission has a flux that is nearly constant with frequency (/3 « —2.1), 
so free-free emission becomes relatively more important than synchrotron at higher frequencies. 
Typically the crossover frequency is near 60 GHz at higher latitudes, but can be 20 GHz or lower 
in specific regions in the Galactic plane. Frequencies above ~60 GHz begin to probe the tail 
(/3 « 2) of vibrational dust emission, which is dominant around 90 GHz. In addition to these three 
foregrounds, much recent work has focused on the possibility of significant emission from rapidly 
rotating dust grains; this emission is thought to peak somewhere in the 10-30 GHz range and fall 
off roughly exponentially at higher frequencies. 

The spectral behavior for diffuse foregrounds is of great interest. The spectrum for synchrotron 
radiation follows the energy distribution for high-energy electrons, which is not a pure power law. 
The highest energy electrons lose energy more quickly and thus are reduced in regions where they 
have not been replenished. Such energy loss shows up as a gradual steepening ( d/3/du < 0) in the 


14 In this paper we use the notation that flux density is S ~ v a and antenna temperature is T ~ i/ 3 , with the 
spectral indices related by (3 — a — 2. Unless otherwise noted, results will be expressed in antenna temperature. For 
the pixel size most commonly used in this work (0?92 x 0?92), the conversion from antenna temperature to flux is 
approximately 4.0(i//22.5 GHz) 2 Jy mK -1 (Page et al. 2003). 
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power law index by about 0.5 at frequencies above 10-100 MHz. Further, while the overall index 
as extrapolated from lower frequencies is (3 « —2.7 (Reich & Reich 1988; Lawson et al. 1987; Reich 
et al. 2004), higher frequencies may preferentially sample more energetic electron populations and 
thus have a flatter index (/? ~ —2.5) (Bennett et al. 2003). Observations of both discrete sources 
(Green 1988; Green & Scheuer 1992) and external galaxies (Hummel et al. 1991) show a wide 
variety of synchrotron behavior. Free-free emission also does not follow a strict power law, but the 
physics is well understood and the variation of the power-law index over WMAP’s bands is so small 
that it can be neglected. Finally, the Rayleigh- Jeans tail for vibrational dust emission (i.e. below 
~100 GHz) has never before been accurately measured and the relevant material properties of the 
dust grains themselves are not fully understood (Agladze et al. 1994; Meny et al. 2007). 

The main focus in this work is on foreground emission. Section 2 describes updates to masks 
and foreground-fitting procedures used in previous WMAP analyses (Bennett et al. 2003; Hinshaw 
et al. 2007). A new method to explicitly marginalize over foregrounds for the low multipole analysis 
is described in a companion paper (Dunkley et al. 2008). A new fitting process is described in 
Section 3, which has the following features: 

• The fitting is entirely in real-space with no spherical harmonic decomposition for any com- 
ponent. 

• The spectral indices of the synchrotron and dust emission are not generally assumed to be 
constant and are allowed to vary down to the scale of the fit (approximately one square 
degree). 

• The fit includes the CMB and automatically generates the full likelihood (including covari- 
ance) for all foreground parameters. 

• The polarization data are included and fit simultaneously with the total intensity data. 

This is similar to the technique of Eriksen et al. (2007), however we fit the CMB in pixel-space, use 
less smoothing on the maps, and attempt to obtain more information about individual foregrounds. 

The results of the fit are described in Section 4. While the fitting technique used here delivers a 
CMB map with error-bars, the map itself has not proved to be any better for cosmological analysis 
and so far has been used only as a check. Implications of the fit are discussed in Section 5. WMAP ’ s 
cosmological results do not depend on the fitting process used here. 


2. Five-year Foreground Fits 
2.1. Masks 

The diffuse foreground masks are updated for the five-year data analysis. The primary reason 
is to mask out free-free emission in the areas of the Gum Nebula and p Oph, while keeping a simple 
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method that applies to the whole sky rather than being ad hoc for these regions. 

The new masks are based on three-year public WMAP data products, specifically the three- 
year K and Q band- average maps smoothed to one-degree resolution. These maps are converted 
to foreground-only maps by subtracting the three-year Internal Linear Combination (ILC) map. A 
cumulative histogram is made of the pixels in each foreground map, which serves as a lookup table 
to find a flux level used to define a cut over the desired percentage of the sky. 

Cuts are made at intervals of 5% in the proportion of sky admitted by the resulting mask. 
The K and Q band cuts at each percentage level are combined. Resulting masks are inspected and 
compared with the masks used in the one and three-year WMAP data analyses. We replace the 
old Kp2 mask with the combined K and Q 85% masks. This is the nominal mask for temperature 
data analysis and is denoted KQ85. We replace the old KpO mask with the combined K and Q 
75% masks (KQ75). 

Each of the chosen masks is further processed by omitting any masked “islands” containing 
fewer than 500 pixels at HEALPix (Gorski et al. 2005) N s ide of 512. Each mask is then combined 
with a point source mask, which has been updated from that described in Bennett et al. (2003) 
and Hinshaw et al. (2007) to include 32 newly detected sources from a preliminary version of the 
WMAP five-year point source catalog. Six sources in the final five-year catalog are not included; 
these are relatively weak, with fluxes of 1 Jy or lower in all WMAP bands. The last step combines 
each mask with the five-year processing cut used to omit the Galactic plane from the mapmaking. 

The three-year polarization mask was based on a cut in K-band polarized intensity combined 
with a model of the dust component (Page et al. 2007). The five-year polarization analysis mask 
is the same as the three-year version, with the exception that it is combined with the five-year 
processing cut. 

The MCMC fit described below uses a version of the combined K and Q 95% mask (denoted 
KQ95, and which is similar to the old Kpl2 mask) to distinguish “inside” from “outside” the 
Galactic plane. The mask was enlarged to account for smoothing, leaving approximately 91% of 
the sky. 


2.2. Internal Linear Combination Method 

The Internal Linear Combination (ILC) method is used to produce a CMB map that is in- 
dependent of both external data and assumptions about foreground emission. By construction, it 
leaves unchanged the component that has the spectrum of the CMB and acts as a foreground fit by 
filtering out the combined spectral shape that causes the most variance in the data. As a minimum 
variance method the ILC is guaranteed to produce a map with good statistical properties, but the 
level of remaining contamination can be difficult to assess. 

The algorithm used to compute the WMAP five-year Internal Linear Combination map is the 
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same as that described in the three-year analysis (Hinshaw et al. 2007). We retain the same number 
of regional subdivisions of the sky and their spatial boundaries remain unchanged from the previous 
definitions. The frequency weights for each region are somewhat different, however, reflecting the 
five-year updates to the calibration and beams. The new ILC regional coefficients are presented 
in Table 1, and the map itself is available on the LAMBDA web site 15 . The coefficients describe 
a filter that nulls certain spectral shapes. A slice in parameter space of the spectra nulled by the 
ILC is shown in Figure 1. 


2.3. Maximum Entropy Method 

The maximum entropy method (MEM) is a spatial and spectral fit using templates that are 
intended to distinguish different low-frequency emission sources. By design, the MEM reverts to 
templates made from external data where WMAP’s signal is low. One of the main goals for the 
MEM was to use high-signal regions to investigate the spectral properties of the foregrounds. The 
error properties for MEM maps are complicated and the model is essentially under-constrained so 
there is no meaningful goodness-of-fit statistic. The MEM maps were not used for analysis of the 
CMB itself. 

The five-year MEM analysis is largely unchanged from the three-year analysis (Hinshaw et al. 
2007). As before, the analysis is done on sky maps smoothed to a common resolution of 1° full width 
at half maximum in all bands. To improve the signal-to-noise ratio, we now use maps degraded to 
HEALPix iV side = 128 pixelization instead of AT side = 256 (the pixel size for the former is 0.46°). In 
the first year and three-year analyses, the logarithmic term that forces the solution to converge to 
the priors for low S/N pixels was missing a factor of e (Cornwell et al. 1999); this has been fixed. 
The model is fit for each pixel p by minimizing the functional H = A + A B (Press et al. 1992), 
where A is the standard x 2 of the model fit, and we now use B = T c (p) ln[ T c (p)/(eP c (p))]. Here 
T c (p) is the model brightness of emission component c (synchrotron, free-free, dust) in pixel p, and 
P c (p) is the prior estimate of T c (p). The parameter A controls the relative weight of A (the data) 
and B (the prior information) in the fit. An iterative procedure is followed that uses residuals from 
the fit at each iteration to adjust the spectrum of the synchrotron component for each pixel. Hence 
any “anomalous” component, such as electric dipole emission from spinning dust, will be included 
in the synchrotron component. See Hinshaw et al. (2007) for a complete description of the analysis, 
including the form of the model, the adopted priors, and the iterative procedure. The MEM maps 
are available for public download on the LAMBDA web site. Figure 2 shows the spectrum of the 
MEM components compared to the total observed foreground spectrum for 20° < |6| < 30°. 

Comparison of MEM results from the five-year and three-year analyses shows an increase in 
the model brightness of all foreground components at high Galactic latitudes. The changes are 


15 http : //lambda . gsf c . nasa . gov/ 
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Fig. 1. — A slice in parameter space of the surface nulled by the ILC coefficients, assuming a three- 
component foreground model with power-law spectral behavior, T{u) = T s u^ 3 + TfU 2Ai + T d v^ d . 
Each line is for a single ILC region, denoted by number. The parameter space is Tf/T s , T d /T s , (3 S , 
0 d . For this plot the x-axis is 0 S and the y-axis is T d /T s . The parameters Tf/T s and 0 d are fixed 
at 0.7 and 1.8, respectively. Each color is a different ILC region. Despite the variety amongst ILC 
coefficients, they often null similar regions of parameter space. 
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mostly due to differences in the zero levels of the five-year and three-year maps. The inclusion of 
the factor of e in the MEM functional also leads to a small contribution. The method of setting 
map zero levels has not changed since the first year analysis. The internal linear combination CMB 
map is subtracted from the 1° smoothed map in each frequency band, and the zero level is set such 
that a fit to the residual map of the form T(|6|) = T p esc |6| ± c, over the range —90° < b < — 15°, 
yields c = 0 (Bennett et al. 2003). The three-year analysis procedure was done using a preliminary 
three-year ILC map in which the monopole was nonzero. Offsets of 21.1, 19.4, 19.3, 19.4, and 19.6 
fiK should be added to the three-year K, Ka, Q, V, and W band maps, respectively, to give maps 
that yield esc |6| fit intercepts of zero when the final three-year ILC map is subtracted. 

Available foreground templates are expected to trace the distribution of foreground emission 
more reliably than a esc |6| model, so template fitting has been done to check the zero levels of the 
five-year maps. The five-year ILC map was subtracted from the five-year 1° smoothed maps, and 
the residual map for each band was fit to a linear combination of synchrotron, free-free, and dust 
templates plus a constant offset. Uncertainties in the zero levels of the templates were propagated 
to obtain an uncertainty in the derived offset value. For the synchrotron template, the 408 MHz 
map of Haslam et al. (1982) was used with an offset of 5.9 K subtracted (Lawson et al. 1987). The 
quoted zero level uncertainty of this map is ±3 K (Haslam et al. 1982). For the free-free template, 
the composite all-sky Ha map of Finkbeiner (2003) was used, with a correction for extinction 
(using the dust extinction map of Schlegel et al. 1998) assuming the dust is coextensive with the 
emitting gas along each line of sight (Bennett et al. 2003). The adopted zero level uncertainty is 
±1 Rayleigh, as estimated by Finkbeiner for the southern Ha data. For the dust template, the 94 
GHz emission predicted by model 8 of Finkbeiner et al. (1999) was used. The adopted zero level 
uncertainty is ±0.2/iK, propagated from a zero level uncertainty of ±0.044 MJy sr _1 for the 100 
fim dust map of Schlegel et al. (1998). 

Fits were done to N s i^ e = 512 pixels that are outside of the combined KQ85 plus point source 
mask and have optical depth at Ha less than 0.5, based on the Schlegel et al. (1998) extinction 
map. This pixel selection covers 74% of the sky. The offsets from the fits are —25 ± 19, —5.4 ± 6.8, 
—2.2 ± 3.9, —2.2 ± 1.5, and —1.5 ± 0.7/zK in K, Ka, Q, V, and W bands, respectively. Thus there 
is no evidence for significant error in the five-year map zero levels as determined from the esc |6| 
fitting. For comparison, northern hemisphere esc |6| fits can be used to estimate uncertainties in 
the zero levels; the northern hemisphere gives offsets of —9.2, 3.2, 3.5, —2.5, —5.9 fiK for K, Ka, Q, 
V, and W bands, respectively, relative to the zero levels calculated from the southern hemisphere. 


2.4. Template Cleaning 

The foreground template subtraction technique used in the five-year analysis is unchanged 
from that used in the three-year release. The method is described in Hinshaw et al. (2007) for 
temperature cleaning and Page et al. (2007) for polarization cleaning; details are not repeated 
here. 
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Fig. 2. — Comparison of the observed foreground emission spectrum (diamonds) with spectra 
of the total MEM model and the individual model components (line segments between WMAP 
frequencies), averaged over 20° < |6| < 30°. The mean model brightness exceeds the mean observed 
brightness at the higher frequencies because the observed brightness is negative for some pixels and 
the model is constrained to be positive for each pixel. This is less apparent in the five-year results 
than in the three-year results because with the higher signal- to- noise ratio and higher map zero 
levels there are fewer negative pixels. 
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In summary, for temperature cleaning a model of the foreground emission is computed from 
a simultaneous fit to the eight five-year Q, V and W-band maps, and that model is then used to 
produce foreground-reduced maps suitable for cosmological studies. 

The model takes the form 

= b\{y)(TK{p) ~ Txaip)) + M^Khc*(p) + h(u)M dust (p) (1) 

where p indicates the pixel, the frequency dependence is entirely contained in the coefficients bi , and 
the spatial templates are the WMAP K-Ka temperature map, the Finkbeiner (2003) composite Ha 
map with an extinction correction applied, and the Finkbeiner et al. (1999) dust model evaluated 
at 94 GHz. All of these spatial templates are available on LAMBDA. 

The coefficients of the model fit to the five-year data are presented in Table 2. Small changes in 
the five-year coefficients compared to the three-year values (Table 5 of Hinshaw et al. 2007) reflect 
the five-year updates to absolute calibration and beam profiles. 

For polarization cleaning the maps are degraded to low resolution (N s - lde = 16). The model 
has the form 

[Q(y,P),U(y,p)] model = ai(v)[Q(j3),U(p)] K + a 2 (v)[Q(p),U(p)} dus t (2) 

The templates used are the WMAP K-band polarization for synchrotron, and a low resolution 
version of the dust template used above with polarization direction derived from starlight measure- 
ments. The coefficients of the model fit to the five-year data are in Table 3. As was the case for 
the three-year data, a fit fixing the synchrotron spectral index was found to have no influence on 
cosmological conclusions and was not used for analysis. 


3. Markov Chain Monte Carlo Fitting 
3.1. Description 

The analysis is carried out with band-averaged maps at each frequency, which are calibrated in 
antenna temperature, smoothed to a one-degree Gaussian beam, and pixelized using an N s [ de = 64 
HEALPix grid. This makes the fit computationally manageable and ensures that pixel-pixel corre- 
lations are small, simplifying the error description. The maps use the esc |6| fit process described 
above to determine the zero-point. 

Next we parameterize the emission in each pixel with a physical model. The model depends 
on the parameters in a non-linear way and the parameters can be highly correlated. A Monte 
Carlo chain is run for each pixel to determine the probability distribution for the parameters of the 
model using the Markov chain technique (Gilks et al. 1996). Because of parameter correlations, 
the matrix describing the optimal step size is not diagonal. The starting points and initial step 
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Table 1. ILC coefficients per region a 


Region 

K-band 

Ka-band 

Q-band 

V-band 

W-band 

0 

0.1336 

-0.6457 

-0.3768 

2.2940 

-0.4051 

1 

-0.0610 

-0.1327 

-0.1873 

1.7691 

-0.3880 

2 

0.0037 

-0.2432 

-0.3792 

1.7956 

-0.1768 

3 

-0.1104 

0.2395 

-0.6424 

1.5032 

0.0101 

4 

-0.0843 

0.1271 

-0.4584 

0.9739 

0.4417 

5 

0.1918 

-0.7238 

-0.4902 

2.4844 

-0.4622 

6 

-0.1052 

0.2614 

-0.6223 

1.0253 

0.4407 

7 

0.0913 

-0.3849 

-0.6033 

2.3288 

-0.4319 

8 

0.2208 

-0.5436 

-1.0938 

3.2084 

-0.7918 

9 

-0.0922 

-0.0695 

-0.1810 

1.2619 

0.0808 

10 

0.1724 

-0.9608 

0.0350 

2.6456 

-0.8923 

11 

0.2374 

-0.8975 

-0.4897 

2.7246 

-0.5747 


a The ILC temperature (in thermodynamic units) at pixel p 
of region n is T n (p) = Yli=i C n,iT l (p), where £ are the coeffi- 
cients above and the sum is over WMAP’s frequency bands. 


Table 2. Template cleaning temperature coefficients 


DA 

b 1 

b 2 (/zK R" 1 ) 

b 3 


hf f h (pK R- 1 ) 

Q1 

0.245 

0.981 

0.201 

-3.18 

5.99 

Q2 

0.243 

1.009 

0.199 

-3.22 

6.01 

VI 

0.058 

0.666 

0.461 

-3.44 

6.38 

V2 

0.056 

0.647 

0.477 

-3.43 

6.38 

W1 

0.000 

0.398 

1.262 


6.62 

W2 

0.000 

0.393 

1.277 


6.62 

W3 

0.001 

0.398 

1.242 


6.61 

W4 

0.000 

0.395 

1.271 


6.62 


a Derived from b\\ W-band values are less than -4. 

b Free-free to Ha ratio at K-band, as derived from b\ and 
62 . The expected value for an electron temperature of 8000 K 
is 11.4 pK R -1 (Bennett et al. 2003). 
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proposal matrices are generated using a “best guess” from the data. In cases where the initial 
guess turns out to be poor, the fitting process is retried using the existing chain to improve the 
guess. Any retries or poorly conditioned proposal matrices are flagged. Each chain is checked for 
convergence using the criteria described in Dunkley et al. (2005), and any lack of convergence is 
also flagged. 

The basic form of the model for each pixel is 


/ V / V \ 0f ( V \ 0d 

T « = T -U) +t '{- k ) 

for the antenna temperature and 

/ v \ AM / u \ A 

Qiy) = Qs ( I “i - Qd ( ) "i - O'iy^Qcmb 

\VK J \VW J 

( V \ 0s(y) ( V \ 0d 

u(u) = u a ( — ) +U d ( — ) + a{u)U cmb 

\VK J \vw J 


(3) 

(4) 

(5) 


for Stokes Q and U parameters. The subscripts s,/, d stand for synchrotron, free- free, and dust 
emission, vk and vw are the effective frequencies for K and W bands (22.5 and 93.5 GHz), and 
a(y) accounts for the slight frequency dependence of a 2.725 K blackbody using the thermodynamic 
to antenna temperature conversion factors found in Bennett et al. (2003). 

For each pixel, the \ 2 °f the fit is then calculated in the standard way 


X 


2 


E D - N - lD ' 


( 6 ) 


where is the difference between the data vector (T, Q, U) and the model vector at each 
frequency. The matrix N*, is the noise covariance matrix, and is derived directly from the N 0 \> s 
maps (rebinned from an N s \^ e of 512 to 64), with minor modifications discussed below. 

Not all parameters in the model are free to vary; a(y) and /3f are fixed by known physics, the 
Q and U parameters for foregrounds are related using K band as a template for the polarization 
angle, and for most of the following f3 s (u) is assumed to be constant with frequency (though allowed 
to vary spatially). The free-free index was fixed at (3f = —2.14; typical variation in this value at 
WMAP frequencies is ±0.015 (Oster 1961, see also Quireza et al. 2006; It oh et al. 2000 for recent 
refinements) which is too small for WMAP to detect. Similarly, results of the fit were not found 
to depend strongly on whether U s and Ud were treated as freely independent parameters. 

The “base” fit which allows for spatially varying synchrotron and dust spectral indices has 
10 independent parameters per pixel: T s , T/, T d , T cmb , /3 a , (3 d , Q s , Q d , Q cm b , and U crnh . More 
restricted fits fixed (3 S , / 3 d , or both. Other fits allowed for a frequency-dependent (3 S by defining 


Psiy) = Ps+c ln(u/u K ) 


(7) 
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where the new parameter c can introduce a gradual steepening (or shallowing). 

Note that the models used here assume that polarized and unpolarized synchrotron emission 
have the same spectral behavior. While this assumption appears to be safe at high latitudes, it may 
not be accurate for lines of sight that pass through the Galactic plane. This is further explored in 
Section 5. 


Finally, some fits allowed for an additional independent component, either using the exact 
“cold neutral medium” (CNM) spectrum for spinning dust (Draine & Lazarian 1998), or using a 
generalized analytic form 


Tsdiy) — A s d 


H ^) 0d+1 

ex P(*V v sd ) ~ 1 


( 8 ) 


The analytic form is a modified blackbody with the amplitude, low-frequency spectral index, and 
turnover frequency explicitly decoupled from one another. Plots of both the exact and “shifted” 
spectra used in the fitting process are shown in Figure 3, as well as a curve showing that the 
analytic form is indeed a good approximation to the numerically calculated spectrum. In practice 
the low-frequency spectral index is irrelevant because the desired shape for foreground fitting has 
u s d well below 22 GHz and is thus dominated by the exponential cutoff. While this form was 
motivated as an analytic approximation to spinning dust spectra it could also represent a variety of 
other physical sources of microwave emission. This component was assumed to have no significant 
polarization. 


Each chain itself is a multi-step process. The code makes an initial guess for the best-fit 
parameters and runs for a burn-in period to find the region of parameter space near minimum y 2 . 
There is then a “pre-chain” to find the approximate moments of the likelihood; these moments are 
used to optimize the proposal distribution for the final chain. Problems at any stage due to lack of 
convergence or poorly characterized parameter distributions are flagged and recorded; only rarely 
are more than 0.5% of pixels so affected, and most problems are due to random fluctuations and 
disappear with longer chains. 


3.2. Tests and Sources of Error 

The Monte Carlo process (with Metropolis steps) has the advantage that it can sample the 
full parameter space and will converge on the likelihood even if the likelihood is non-Gaussian 
or unknown a priori. The disadvantage is that degeneracies in parameter space will slow the 
convergence, and cutting off regions of parameter space to improve convergence can bias the results. 
The prime example of this is degeneracy between synchrotron and free-free emission amplitude. If 
the synchrotron spectral index is allowed to flatten to the free-free value then the amplitudes of 
the two components become degenerate parameters, which can distort the fit. 

We test this with simulated maps where the input foreground is known. To ensure that the 
noise properties of the maps are well understood we used extensive simulations. We combined the 
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Fig. 3. — The exact cold neutral medium (CNM) spinning dust spectrum as calculated by Draine 
& Lazarian (1998), an analytic fit to the model, and the ad hoc “shifted” model which better fits 
radio observations in the Galactic plane. The shifted model may represent a mix of warm neutral 
medium and warm ionized medium models, or another emission process entirely. The vertical axis 
is in units of antenna temperature, but the overall scale is arbitrary. Agreement between the exact 
model and the analytic approximation is better than 5% over the frequency range from 15 to 35 
GHz. This is smaller than the fractional error from fits that include a spinning dust component. 
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high resolution noise information with the low resolution pixel-pixel covariance in order to generate 
noise realizations that are as realistic as possible, which are then smoothed using the same process 
used for the real sky maps. We then produce mock sky maps with CMB realizations synthesized 
from the WMAP “best-fit” ACDM model, noise realizations from the five year noise covariance 
matrix, and foregrounds generated from a variety of models. 

When degeneracies exist the random fitting process tends to share the amplitude evenly be- 
tween degenerate parameters. This can lead to biasing if the true sky does not also have equal 
contributions from such parameters. This effect can be seen clearly by comparing the reconstruc- 
tion of the synchrotron and free- free components (Figure 4). Figure 5 shows histograms of a single 
pixel chain from the mock sky fit, with the input values of the parameters marked with a bold 
cross. 

The fitting process uses the error information contained in the iV 0 bs maps which includes 
covariance between the Q and U Stokes parameters within a pixel, but covariance between pixels 
(due to low-frequency noise or the smoothing process) is not included. The correlation coefficient 
between adjacent pixels ranges from less than 0.20 up to 0.45 (for K-band and W-band, respectively) 
due to the smoothing process. The fit treats pixels as independent, essentially marginalizing over 
all other pixels when fitting, so the main effect of the correlations is to introduce similarly small 
pixel-pixel correlations in the y 2 values. This has a negligible effect on the results as long as 
goodness-of-fit is averaged over large enough regions. 

Smoothing also reduces the overall noise level, and this has been modeled through direct 
noise simulations and accounted for in the fit process. The method used was to generate many 
realizations of simulated noise maps based on the original iV 0 b s information, smooth them using 
the same window functions used to smooth the real data, and fit a Gaussian shape to a histogram 
of the result. From this process an overall multiplicative rescaling factor was determined for each 
frequency band and applied to the iV 0 bs files used for the final fit. With this correction, the y 2 per 
degree of freedom should be close to unity for an ideal fit; for the mock fit described above the 
mean reduced y 2 was 1.11. 

Small differences in the beam solid angle from one frequency band to another can distort the 
inferred spectral index, especially near bright sources. We used the Jupiter-based beam maps from 
Hill et al. (2008) and smoothed them to a common one-degree beam, similar to how the sky maps 
are smoothed, and found that beam systematics at the known level can cause deviations of up to 
±0.1 in the spectral index. Errors of this type (that are multiplied by the sky signal) are included in 
the fit by adding 0.3% of the antenna temperature to the error budget for each pixel. This number 
was derived from the observed scatter one beam- width away from bright point sources. Systematics 
of this type do not average down, and so can quickly become dominant at low resolutions. 

The smoothing kernel used to match the bands to one-degree resolution uses the symmetrized 
beam profile, and hence does not take into account beam asymmetries. WMAP ' s observational 
strategy, however, symmetrizes the beam to a large extent. Page et al. (2003) investigated the 
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Output fit - input mock sky (mK) 


Fig. 4. — Difference maps between input and output foregrounds for the mock sky reconstruction. 
For comparison, the peak temperature for synchrotron plus free- free inputs was 17 mK, and for 
the input dust the peak was 1.9 mK. The degeneracy between synchrotron and free- free emission 
means their sum is much better constrained than either component individually. The scatter is 
larger than the bias, which is small compared to the input signal and within the error estimate 
given by the fit. These degeneracy issues are also illustrated in Figure 5. Synchrotron and free- free 
antenna temperatures are defined as measured at K-band, and dust as measured at W-band. 
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Fig. 5. — Histograms for foreground parameters for a single pixel of mock sky reconstruction. 
Dashed lines indicate the mean of the chain, crosses indicate the best-fit point, and stars indicate 
the input values for the parameters. The strong T s -Tf degeneracy and the curved T S -/3 S degeneracy 
are typical for high signal-to-noise pixels. Other parameter correlations are minor. Both the mean 
and best-fit parameters of the chain are within the expected error range from the input values. 
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extent to which remaining beam asymmetry could affect the beam window functions and found it 
to be < 1%. Any effect on the maps due to beam asymmetries should be weak near the ecliptic 
poles and for extended emission not aligned along the plane of the ecliptic. 

There remain small uncertainties of a few (iK both in the true zero-level of the maps and in 
the dipole subtraction process. The offsets primarily affect foreground estimation by changing the 
apparent spectral index when averaging over large, very low signal-to-noise regions. We avoid this 
by explicitly de-weighting or masking pixels with weak spectral index constraints when reporting 
results. Second, as a purely pixel-based method, the MCMC foreground fitting process we use is 
free to produce foreground (or even CMB) maps with non-zero monopole and dipole contributions. 
The sensitivity of the fit to offsets was checked by adding offsets of 100 iiK (several times larger 
than the error as estimated from the esc |6| fits) to the sky maps and repeating the analysis. No 
foreground component was found to change by more than 10% for pixels where the signal-to-noise 
was significant. 


4. Fit Results and Comparisons 
4.1. MCMC Fit 

Each pixel fit consists of 15 data points (Stokes I, Q, and U for each of the five frequency 
bands) and a foreground fitting model can use from 8 to 12 parameters per pixel. The fitting 
process produces a x 2 value for each pixel. Normally the reduced x 2 is found by dividing by the 
number of degrees of freedom. However, the true number of degrees of freedom in this case is 
difficult to determine because neither the data points nor the fitting parameters are statistically 
independent of one another. Using the MCMC chain for each pixel, though, it is possible to use 
the “Bayesian complexity” (described in a cosmological context in Kunz et al. 2006), defined as the 
difference between the average x 2 over the chain and the x 2 of the best fit. This serves as a measure 
of the effective number of degrees of freedom, and can then be used to determine the reduced x 2 per 
pixel. Using the simulated skymaps described above, we have found that the statistical behavior 
of the reduced x 2 defined this way is consistent with that of a x 2 distribution. This “effective” 
reduced x 2 is how we quantify the goodness-of-fit in the tables and figures. 

Pixels with high reduced x 2 are not being well-fit by the model. Since such pixels are largely 
confined to the plane, the sky was divided into regions “outside” and “inside” the Galactic plane 
by using progressively smaller masks until the average “outside” x 2 w as no longer independent 
of the mask. Regions near known point sources from Wright et al. (2008) are excluded from all 
analysis, both inside and outside of the plane, leaving 92% of the full sky. Flagged pixels are 
also not included; for most fits such pixels arise from poorly conditioned covariance matrices, are 
uniformly distributed, and make up less than 0.5% of the sky. 

The x 2 results of several fits are shown in Table 4. Histograms for a single pixel chain are 
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Table 3. Template cleaning polarization coefficients 


Band 

ai a 

Paiy K ,v) h 

«2 a 

Pd(v, vw) h 

Ka 

0.3161 

-3.17 

0.0165 

1.35 

Q 

0.1765 

-3.04 

0.0147 

1.85 

V 

0.0595 

-2.96 

0.0366 

1.58 

W 

0.0450 

-2.35 

0.0822 



a The di coefficients are dimensionless and pro- 
duce model maps from templates. 

b The spectral indices refer to antenna tempera- 
ture. 


Table 4. Model fits to WMAP temperature and polarization data 



# of 

Best-fit x 2 a 


model 

params 

outside plane b 

inside plane b 

full sky 

base 

10 

1.14 

2.23 

1.24 

base + Haslam 

10 

1.14 

2.36 

1.26 

loose priors 

8 

1,09 

3.26 

1.29 

steep 

10 

1.14 

0.97 

1.13 

exact sd 

9 

1.21 

1.63 

1.25 

shifted sd 

9 

1.24 

1.00 

1.22 

ft = -3-2, ft = 1.7 

8 

1.16 

4.33 

1.45 

ft = -2-6, fi d = 1.7 

8 

1.30 

3.42 

1.50 

fi s variable, fid = 1-7 

9 

1.16 

2.92 

1.32 

fi s variable, fid variable 

10 

1.14 

2.23 

1.24 


a Reduced x 2 averaged over pixels in the region, with effective degrees of 
freedom determined by the chain. The statistical errors are less than 0.01. 

b The mask used to define these regions is a smoothed version of the 95% 
mask, the 5- year release analogue of the Kpl2 mask. 
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Fig. 6. — Histograms for foreground parameters for a single pixel of the observed sky from actual 
WMAP data, not simulation. Dashed lines indicate the mean of the chain, and crosses indicate 
the best-fit point. The observed sky shows the same basic behavior as the simulated sky used for 
testing. 
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Fig. 7. — Temperature maps for foreground components as determined by the MCMC fitting 
process for the “base” model. Maps from other models are qualitatively similar. Synchrotron and 
free- free temperatures are as measured at K-band, dust is measured at W-band. Gray pixels are 
those masked due to point sources or flagged as problematic. Top - synchrotron; middle - free-free; 
bottom - dust 
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Fig. 8. — Residuals of the “base” model foreground fit subtracted from the data. The left column 
shows the residuals in units of noise sigma, from -6 to 6. The right column has the 85% mask 
applied and a scale of -3 to 3. Frequency bands are (from top to bottom) K, Ka, Q, V, W. The 
main feature is that the model underestimates Galactic flux in Ka band and overestimates it in Q 
band by a factor of 3 to 5 times the pixel noise. Outside the analysis mask the residuals to the 
model are consistent with noise at the expected level. 
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Fig. 9. — Difference between the MCMC fit and MEM maps. The top panel is the difference (at K 
band) in synchrotron plus free-free emission. Differences between these two components separately 
are larger due to the degeneracy direction. The bottom panel is the difference (at W band) in the 
dust maps. The differences are roughly one percent of total emission at K band and a few percent 
at W band. 
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shown in Figure 6. The mock sky simulations appear to capture the basic behavior of the parameter 
correlations. Detailed comparisons of specific results are in the subsections below. Figure 7 shows 
the basic results and Figure 8 shows the temperature residuals of the best-fit base model subtracted 
from the data, in units of one-sigma of noise. Figure 9 shows the difference between the MCMC 
fit and the five-year MEM fit. The overall fit outside the complex and troublesome Galactic plane 
region approaches a xt of 1*14 and the residuals are mostly randomly distributed, which suggests 
that the overall fit works reasonably well and that the noise properties have been described properly. 

The “base” model uses the 10 parameters described above. Another fit is done including data 
from 408 MHz (Haslam et al. 1981), assuming 10% calibration errors. As a check, a “loose priors” 
fit is done which allows foreground temperatures to become negative. For this fit to converge the 
spectral indices must be fixed, however, so it only uses 8 parameters. The “steep” model fixes 
the dust spectral index but allows for a synchrotron steepening parameter c as described above. 
The “exact sd” model uses the CNM spinning dust spectrum for an additional foreground, whereas 
the “shifted sd” model uses the generalized spectrum with — 4.9 GHz. In both cases, the 
synchrotron and dust spectral indices are fixed to be the same at each pixel (again for convergence 
reasons). Finally, several fits were done with fixed spectral indices to examine the effect of using 
different values, shown as the last part of Table 4. The last entry of the table is the “base” model, 
repeated for ease of comparison. 


4.2. Overall Foreground Features 

Figures 7-12 show maps of the results from the “base” fit. Figures 7 and 10 show the three 
foregrounds themselves and their errorbars. The maps are in units of antenna temperature as 
measured at K-band for synchrotron and free-free emission, and at W-band for dust emission. 
Figure 11 shows spectral index maps binned to lower resolution, where color indicates the spectral 
index and the size of the circle indicates the significance of the fit result at that location. Figure 
12 shows the best xt value achieved at each pixel on the sky. 

Almost regardless of foreground model the fit works quite well outside the Galactic plane, 
giving low xl and foreground maps that are in good agreement with the MEM templates (Hinshaw 
et al. 2007) and other works (Eriksen et al. 2007). Error maps for synchrotron and free-free emission 
have similar morphology due to the degeneracy between their amplitudes. 

The overall dust brightness seems to be largely consistent with the template prediction (Finkbeiner 
et al. 1999), though the fit appears to prefer a spatial distribution somewhat less sharply peaked 
toward the Galactic center (Figure 13). The excess of observed emission compared to that predicted 
at 90 GHz, seen in the original model comparison with COBE DMR data, is still present. Since 
the fit in the plane has high x 2 and is untrustworthy, the overall preferred spectral index for dust 
may be < 2.0, but the significance of this is not high. Weighting with the covariance information 
from the fit and masking low-signal regions, the average value for (3d in the “base” fit is 1.8 with 
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Fig. 10. — Error maps for foreground components as determined by the MCMC fitting process. 
Error maps for synchrotron and free-free emission are similar due to the parameter degeneracy 
between them. Synchrotron and free-free temperatures are as measured at K-band, dust is measured 
at W-band. Gray pixels are those masked due to point sources or flagged as problematic. Top - 
synchrotron; middle - free-free; bottom - dust 
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±0.3 from statistical errors and ±0.2 from systematic error depending on how the cuts are defined. 
For comparison, model 8 of Finkbeiner et al. (1999) predicts fa — 1.55 ±0.01 for a comparable sky 
cut. 


The free-free component is consistent with expectations from previous fits and with Ha ob- 
servations when dust obscuration is taken into account. Free- free emission is quite high in the 
Galactic plane and in several regions (including Gum and Orion) appears to be dominant over 
synchrotron, even in K-band. The ratio of the “base” fit free-free map to extinction-corrected Ha 
map of Finkbeiner (2003) was used to make a map of hff (the temperature-Ha intensity conver- 
sion factor). A histogram was then made of all pixels with intensities larger than 5 Rayleighs (to 
mask out low-signal regions) and less than one magnitude of extinction (using the reddening map 
of Schlegel et al. 1998). A gaussian fit to the peak of the histogram gives hff = 11.8 ± 8.8 (. iK R -1 
at K-band, comparable to the value of 11.4 fiK R -1 expected from an electron temperature of 8000 
K (Bennett et al. 2003), but also consistent with the lower values from template cleaning. 

While subject to degeneracy with the free-free emission, synchrotron radiation is a stronger 
signal in WMAP data than dust emission. Pixel-by-pixel constraints become poor far from the 
plane, however there are still constraints on the best-fit spectral index. For example, by comparing 
fits with constant spectral index, the Northern Polar Spur and the Fan region prefer an index 
of —3.0 or steeper. All fits including total intensity data show the same preference for shallower 
spectral index in the plane as concluded by Bennett et al. (2003); Figure 14 shows fa as a function 
of latitude for a number of different fits. 

From the polarized data, the synchrotron polarization fraction indicates strong depolarization 
toward the Galactic plane consistent with Kogut et al. (2007). Since Faraday rotation should not 
be large at these frequencies, this effect is due to multiple magnetic field orientations along the 
line of sight. Dust polarization fraction appears to follow a pattern similar to the synchrotron 
polarization fraction, though the signal-to-noise ratio is low. This is physically reasonable, as the 
polarization fraction is largely affected by the coherence of the magnetic field along the line of sight. 
This implies that the dust intensity times a constant fraction may not be the best template to use 
for dust polarization in the Galactic plane. 


4.3. The Galactic Plane 

Regions at very low latitudes are not as well fit by the “base” model, and there is dependence 
both on foreground model and fit parameters. A map of poorly fit regions reveals that they are in 
the brightest parts of the Galaxy, where at these frequencies the free-free emission dominates. 

Pixels poorly fit by the “base” model have some common characteristics. Most are bright, 
but this is probably because similar less bright pixels have lower signal-to-noise and thus lower 
X 2 . Many have a K-Ka temperature spectral index similar to what one would expect from free- 
free emission, but a considerably steeper Ka-Q spectral index. Such pixels tend not to be highly 
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polarized, and what polarized emission exists appears to be consistent with synchrotron emission 
with a typical spectral index of (3 « —3. Data for several such individual pixels are shown in Figure 


15. 


For the published cold neutral medium spinning dust model (the “exact” model of Table 4) , 
the maximum fraction of Ka-band flux attributable to spinning dust is 17% outside of the Galactic 
plane (using the 85% mask). The maximum full-sky fraction of Ka-band flux attributable to 
spinning dust is 20% for this model. However, this model still does not provide a good fit within 
the Galactic plane (xl i n this region is 1.63). 

Allowing the spinning dust spectrum to shift in frequency to obtain a better fit results in a 
Ka-band flux fraction of 14% for spinning dust, roughly independent of sky cut. A map of the 
spinning dust component from this fit and its error is shown in Figure 16. The morphology lies 
somewhat between that of dust and free-free emission, though the details depend on the specifics 
of the model. The Galactic plane is equally well-fit by adding a synchrotron steepening parameter 
c into the fit. The actual value of c is generally not well-constrained, but the average value in the 
plane is —1.8. This very rapid steepening does not appear to be consistent with cosmic ray models 
(Strong et al. 2007), but may have some other physical origin. 

Since the goodness-of-fit outside the plane is not improved by the addition of a spinning dust 
component, and low signal-to-noise regions bias the spinning dust fraction upwards, we regard the 
spinning dust fraction of the fits above as an upper bound to the overall amount of diffuse spinning 
dust emission present. As with previous WMAP fits, this new fitting technique continues to find 
that spinning dust is a subdominant emission process. 


5. Discussion 

5.1. Effect on CMB and Cosmology 

The uncertainties of the fit in the Galactic plane preclude CMB analysis for those regions. 
Fortunately, such regions appear to be tightly confined to the plane inside a very narrow sky cut 
(9% of the total sky) and thus can be excluded without losing much information for cosmological 
analysis. The foreground maps from the MCMC fit are similar to those from the MEM fit and 
other foreground templates, which means CMB polarization maps cleaned using such templates 
will also be similar. 

Outside of a narrow band on the Galactic plane the CMB map produced by the fit is visually 
identical to the ILC map. The difference between the “base” fit CMB map and the ILC map is 
shown in Figure 17. The total variance of this difference map outside of the KQ85 mask used for 
power spectrum analysis is 116 (i K 2 , much lower than the CMB power. The variance between 
the “base” fit CMB map and the “shifted sd” CMB map is 44.1 fi K 2 outside of the KQ85 mask; 
variance from one fit to another is generally even smaller for other combinations. 
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The CMB polarization maps produced by the MCMC fit presented here are noisy and show 
some evidence of synchrotron contamination. Nevertheless, the covariance maps from the fits can 
be used to bound the amount of contamination present, and are available on the LAMBDA website. 
For cosmological analysis a different method is used to marginalize over polarization foregrounds. 
For a full description see Dunkley et al. (2008). 


5.2. More Complicated Models 

All of the models so far fit assume simple forms for foreground emission. These are justified if 
the observed emission is dominated by a few emission mechanisms which simply combine additively 
along the line of sight and have minor spatial variation within the beam. In more complex regions 
of the Galaxy, however, things may not be so simple. 

If two synchrotron regions along the line of sight have their polarization angles oriented nearly 
orthogonally, then the total polarized emission will be sharply reduced. If the two regions have 
different spectral indices then the cancellation in polarization will be maximized at the frequency 
where the individual polarization amplitudes match, causing a dip in the polarization spectrum. 
Thus even with pure synchrotron emission the polarization spectrum can look quite different from 
the temperature spectrum. 

To assess this effect 100,000 Monte Carlo realizations were made of a superposition of two 
independent randomly oriented synchrotron emitting regions. Parameters for the distribution of 
intensity and spectral index were chosen to roughly correspond to observations, but the simulation 
was meant only to provide a generous estimate of how different temperature and polarization 
behavior could be simply due to multiple synchrotron regions along the line of sight. The mean 
spectral index difference was small (—0.051) but the standard deviation was not insignificant (0.12) 
and the distribution was non-gaussian with high kurtosis (4.4). Over 20% of the simulations had 
an absolute slope difference larger than 0.1. Thus we cautiously conclude that while differences 
between temperature and polarization spectral indices at the ~0.1 level could be quite mundane 
in origin, consistent differences of 0.25 or larger are probably not due to chance alignments in 
polarization angle and may be caused by an unpolarized non-thermal temperature component. 

For the fit, free- free emission was modeled as a pure power law based on the assumption 
that the plasma is optically thin. In reality, H II regions can become dense enough to become 
optically thick at frequencies as high as 20 GHz, although such regions are spatially small and do 
not contribute significantly to the observed emission for a beam as large as WMAP’s. Further, 
to obtain rising flux at K-band requires very high emission measure (~10 9 cm~ 6 pc). Even if they 
were somehow significant, such regions can not explain a steepening spectrum past Ka-band. 

Synchrotron self-absorption can also cause a low-frequency turnover, but the physical param- 
eters necessary for the turnover frequency to lie in or near the WMAP range imply conditions 
typically only found in active galactic nuclei or other extreme regions. It may be physically possi- 
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ble for synchrotron radiation from stellar-mass black hole jets or accretion disks to become optically 
thick at WMAP frequencies, but this phenomenon has yet to be clearly observed and it is unlikely 
that such emission would contribute significantly at WMAP ' s resolution. 

Synchrotron radiation is suppressed when emitted from a region with a refractive index less 
than unity, such as a plasma. This is known as the Tsytovitch-Razin effect, and causes strong 
suppression of synchrotron emission below 20 (n e /B) Hz, where n e is the electron density (in cm -3 ) 
and B the perpendicular component of the magnetic field (in G). For typical Galactic electron 
densities and magnetic fields, this cutoff is in the 3-300 MHz range, at most. Unless energetic 
electrons play an unexpectedly large role in diffuse emission, WMAP should not see significant 
Tsytovitch-Razin suppression. 

Diffusive synchrotron radiation (DSR) differs from an ideal synchrotron spectrum because of 
the presence of significant random fluctuations in the magnetic field (Fleishman 2005). In this 
model lower-energy electrons experience small-scale turbulence in the magnetic field structure and 
follow non-circular paths due to the random deflections. In such models the emission spectrum can 
turn over from a power law with (3 « —2.1 in the turbulence-dominated diffusive regime to a normal 
synchrotron spectrum at higher frequencies. While in most models this occurs at low frequencies, 
there is some indication that for pulsar wind nebula the turbulence is relevant up to the GHz range 
and above (Fleishman & Bietenholz 2007). Whether DSR can occur for less compact objects is not 
understood at this time. 


5.3. Other Components 

Much has been written on the possible presence of anomalous emission in the lower frequency 
bands of WMAP. There are at least two categories of observations: one is of emission that is diffuse 
and significant over large portions of the sky outside the Galactic plane, and another where the 
emission is important and perhaps even dominant in specific compact regions. 

Much evidence for diffuse anomalous emission comes from template correlations (de Oliveira- 
Costa et al. 1999; Bonaldi et al. 2007; Dobler & Finkbeiner 2007) rather than direct fitting of the 
data (though for a recent example of the latter, see Miville-Deschenes et al. 2008). Characterizing 
the error and offsets in templates made from data at very different frequencies has proven challeng- 
ing. Nevertheless, Dobler & Finkbeiner (2007) show that using an Ha template to fit WMAP data 
results in an improvement of the xl by 0.016 (to 2.977), and that the spectrum has a significant 
“bump”. Though the improvement seems small, due to the large number of degrees of freedom 
it is statistically significant, and appears to be robust against the systematic error investigated in 
that work. The Ka-band excess for select Galactic plane pixels (Figure 15), however, is in the data 
alone, independent of any template. Boughn & Pober (2007) also find from combining 19 GHz 
data with WMAP K-band that the Galactic plane seems to have antenna temperature falling less 
steeply than (3 = — 2 (i.e. a rising flux spectrum). 
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Specifically regarding compact regions, Finkbeiner et al. (2002) previously reported on two 
regions which might show excess emission in the 10-40 GHz range due to spinning dust. CBI ob- 
servations (Dickinson et al. 2006) failed to find anomalous emission from one, LPH96 201.663+1.643; 
other authors (McCullough & Chen 2002) had previously raised the possibility that such emission 
might be due to an optically thick ultracompact H II region. The other, LDN 1622, was found by 
Casassus et al. (2006) to have a spectral energy distribution consistent with spinning dust, a result 
driven in part by the lack of flux at 5 GHz found by the Parkes-MIT-NRAO survey of Condon et al. 
(1993). Other surveys at lower frequencies with a larger angular resolution more comparable to 
WMAP’ s, however, have not measured a lack of flux - the 408 MHz data summarized by Haslam 
et al. (1981) measures more flux than WMAP K-band and is consistent with a mixture of power-law 
thermal and non-thermal components for the region containing LDN 1622. 

Separately, Scaife et al. (2007a) recently observed a sample of northern H II regions and found 
no evidence for anomalous emission in any, but observations with the Very Small Array (Scaife 
et al. 2007b) find some evidence for a 33 GHz excess in SNR 3C396. Thus for compact regions the 
status of anomalous emission appears to be mixed. 

There has been some discussion in the literature of correlation between CMB maps and neutral 
hydrogen (Verschuur 2007), but this result was not found to be statistically significant (Land & 
Slosar 2007). 


5.4. Directions for the Future 

It is quite probable that at least one of the above model complications or additional components 
is relevant for understanding our Galaxy. More data is needed, particularly in the 5 to 30 GHz range. 
Further, the inability to measure flux at large angular scales is a problem for many observations, 
particularly when the angular scale limit depends on the observing frequency. This continues to 
make precise comparison of results difficult. Large-scale observations with calibration errors at 
the percent level or better are needed to address the nature of features seen in some pixels of the 
Galactic plane. 

For the WMAP foreground fits the dust spectrum was treated as a pure power law. In reality, 
dust emission in the WMAP bands is probably dominated by a cold component with a low enough 
temperature that the exponential cutoff is not negligible. However, since the frequency range of 
WMAP ’ s dust sensitivity is narrow the largest effect of the exponential correction is simply a 
modification of the apparent power-law index, which for typical cold dust temperatures (~10 K) 
amounts to a change of about 0.1 in /?^. This effective bias is therefore small compared to typical 
errorbars. However, there are already hints that the extrapolation of dust models to millimeter 
wavelengths is not entirely satisfactory. 
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6. Conclusions 

• WMAP’s temperature and polarization data outside the Galactic plane are well described 
by the standard three foreground components: synchrotron, free-free, and thermal dust, each 
with power-law spectral indices. 

• The spectral index for synchrotron radiation at high latitudes is consistent with (3 ^ — 3, 
with trend toward (3 = —2.7 seen at lower latitudes. The spectral index for dust is not well 
constrained but appears consistent with (3 = 2. 

• Some localized regions in the Galactic plane show emission with (3^—2 below 33 GHz which 
steepens by as much as A (3 = —0.8 above 33 GHz, and this emission is mostly unpolar- 
ized. Both spinning dust and synchrotron steepening models can be used to fit this emission 
component, whose physical origin is unclear. 

• CMB maps from different model fits show < 50 /iK 2 of variance relative to each other outside 
the KQ85 analysis mask, and <120 /iK 2 of variance relative to the five-year ILC map. The 
CMB and cosmological results are robust to changes in the foreground model. 

• WMAP serves as a precise (< 1% error), unbiased, full-sky survey of the Galaxy which can 
reveal large-scale microwave emission features never before seen. 
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NNX07AL75G SOI, LTSA03-000-0090, ATPNNG04GK55G, and ADP03-0000-092. This research 
has made use of NASA’s Astrophysics Data System Bibliographic Services. We acknowledge use of 
the HEALPix, CAMB, and CMBFAST packages. 
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Fig. 11. — Map of synchrotron spectral index for the “base” fit, binned to N s ^ e = 16. Color shows 
the value of the spectral index, and circle area indicates the weight crp ~ 2 given by the fit. Pixels 
with Xu > 2 were explicitly de-weighted. 
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Fig. 12. — Map of the reduced x 2 per pixel achieved by the “base” fit. Except for a slight excess of 
high x 2 values from the plane, the overall statistical distribution is that of an ideal x 2 distribution. 
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Fig. 13. — Residuals between the best-fit dust map from the “base” model and 1.27 times the 
“model 8” dust map of Finkbeiner et al. (1999) extrapolated to 94 GHz. Pixels with > 5 
or covered by the point source mask are colored gray. The MCMC fit prefers less dust emission 
directly toward the Galactic center, and somewhat more emission further away. Units are antenna 
temperature at 94 GHz. 
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Fig. 14. — Synchrotron index plotted as a function of latitude for several fits. Pixels are binned 
by latitude, and only longitudes between 350 and 10 degrees are included. Error regions indicate 
the 68% scatter within each bin. Solid black squares (gray region) are for the “base” model, blue 
crosses (blue region) are for the “base” model with 408 MHz data, and red empty squares (red 
region) are for a fit using only WMAP's polarization data. Fits where the dust spectral index was 
fixed to 1.7 and 2.0 are almost identical to the “base” fit. The trend to flatter spectral index in 
the plane does not appear when only polarization data are used, but the signal-to-noise is not high 
enough for the discrepancy to be significant. 
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Fig. 15. — Temperature (Stokes I) and polarization ( y/Q 2 + U 2 ) spectra for poorly fit pixels in 
different regions of the sky. The vertical scale is arbitrary flux units, but polarization data is shown 
beside the corresponding temperature data with a number in red indicating the approximate (K- 
band) ratio of intensity to polarization. Polarization data is not shown for pixels with low signal 
to noise. Numbers in parentheses are Galactic coordinates. The blue curve is the best-fit “base” 
model. The model does not include synchrotron steepening; convex slopes are due to negative 
CMB contributions. The grid of gray lines indicates spectral indexes of a = —0.14 and a = — 1. 
The left plot shows pixels further from the plane; the right plot shows pixels near the plane. The 
fit converged for all regions shown and none were covered by the point source mask. However, the 
four pixels shown on the left are (from top to bottom) within a few degrees of p Oph, £ Oph, Tau 
A, and the Orion nebula. 
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Fig. 16. — Top: Map of spinning dust antenna temperature T sc l (as measured at K-band) for the 
“shifted” spinning dust model fit. While the model is based on a spinning dust spectrum, the 
data are insufficient for distinguishing the physical emission mechanism from other possibilities 
mentioned in the text. Maps of the other three foregrounds are qualitatively similar to that of 
the “base” model, with the synchrotron amplitude reduced to accommodate the additional low- 
frequency component. Bottom: Map of spinning dust temperature divided by the marginalized 
temperature error as calculated by the MCMC chain. This is statistical error only and does not 
include uncertainties in the model. 
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Fig. 17. — Difference map between the MCMC fit for the CMB and the ILC map. The variance 
outside the KQ85 analysis mask is 116 /iK 2 , lower than the primordial CMB power. 



